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ABSTRACT 


Summary  results  of  U.  S.  Navy  ^Underwater  Sound1  Laboratory  experi¬ 
mental  program  to  investigate  acoustic  phenomena ;in  the  Arctic  Ocean  are 
presented. 

Since  the  sound  velocity  is  an  increasing  function  of  depth,  propaga¬ 
tion  is  characterized  by  upward  refraction  and  surface  reflection  (RSR)., 

A  rough-surface  model  of :,the  ice  coyer  accounts ‘for  both  forward  and 
back-scatter.  The  roughness -wavelength  spectrum  calculated  from  re¬ 
verberation  measurements  is  similar  to  that  for  the  sea  surface,  although 
the  level  is  much  higher.  Forward1  scatter  loss  depends  on,. total  rough¬ 
ness' and  is.-  responsible  for  severe  attenuation  of  high  frequencies.  Propa¬ 
gation  and  reverberation  data  both  imply  an  R.  M.  S.  roughness  of  from 
two  to  three  meters,  which  is  consistent  with  under-ice  profile  measure¬ 
ments.  Propagation  of  explosive;  wavesnis  described  by  normal  mode -and 
ray  theories:  At  short  ranges -convergence  zones  are  observed.  Because 
the  ice  cover  shows  a  "  critical,  angle"  ^dependence,  the  time  disper sion  yf 
the  wave  train  at  long  range  in  deep  water  is -quite  well  defiried  t  =  ’ 

In  shallow  water,  the  bottom  may  produce  bottom  reflection  modes',  y  re¬ 
duced  dispersion  of  the  refractedmode. 

Unusually  low- ambient  noise  levels- are  observed  during  undisturbed 
periods;  Spectra  indicate  that  the  noise  background  arises  mainly  via 
long  range  propagation  at  these  times.  >In  periods  -of,  high ’wiudvahd  parti¬ 
cularly  rapid  temperature  change,  characteristic  ice  noise  save  predomi¬ 
nantly -of  local  origin.  -Occasionally  during  the  summer  -biological  sounds 
of  various  forms  of  marine  life  are  heard.. 
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1.  INTRODUCTION 


1.  1  Brief  History 

In  the  summer  of  1958,  the  Underwater  Sound  Laboratory  ( USL) 
began,  an  experimental  study  of -underwater  sound,  propagation  in  the  Arctic 
Ocean;1*  Signals'  from  underwater  explosives  were  exchanged  between  two 
“drift  stations,  Fletcher's  Ice  Island  (T  3)  and  Station  Alpha,  which  at  that 
time  wer,e  some  SOO^.km -apart.  These  early  experiments  revealed  some 
pfithe  unique  features  of  long  range  propagation  under  the  ice  cover.  Notar 
bly,  -the  arrivals  were  .found  to  consist  of  a  dispersive  quasi-sinusoidal 
wave; train  of  10  —  l'OO-  Hz,  frequency  range.  The  dispersion  and  transmis¬ 
sion  ioss  were  qualitatively,  explained  by  a  half  space  sound  channel'  model 
in  which  the  attenuation- of  high, frequencies  was  ascribed  to  scattering  by 
the  rough  undersurface  of  the  ice. 

During  the  summer  of  1959,  the  experimental  program  was  continued 
between  T  3;  and  Station  Charlie  at  a  range  of  about  1200  km.  2  In  order  to 
study  propagation  as a  function  of  range,  %  series  of  aircraft  flights  were 
made  in -which  practice  depth  charges-,(PDC)  were  dropped  at  intervals. 

The  locations  of  the  stations  and  aircraft  drops  are  shown  in  Figure  1.  Sig¬ 
nals  were  recorded  at' .both  station  sto  provide  data  for  a  variety  of  propa¬ 
gation  paths  of  varying  length,  involving  bothdeep  and  shallow  water.  In. 
September*,  the  visit  of,  the  Icebreaker  Staten-  Island  provided  the  opportu¬ 
nity ‘to  obtain  data  .at  intermediate  ranges  which,  was  difficult., to- obtaimwith 
fast.rhov.lrig  aircraft;  Other  local  measurements, were  made  in  the  vicinity 
of'T  3  out' to  a  range  of-  10  ,km. 

During  April'  and'  May  of  19-62,  a,-co6perative  experiment  by  USL, 
Lament  Geophysical  Laboratory  and  the  Pacific  Naval. Laboratory  was 
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carried  out  between  three  drift  stations:  T  3,  Arlis  II  and  Polar  Pack  I 
at  locations  shown  in  Figure  1.  The  main,purpose  of  this  series  was  to 
obtain  supplementary  cold  weather  data  prior  to  the  summer  warmup. 

Some  aspects  of  the  program  have  been  previously  reported  in  two 
papers,  one  on  propagation*  and  another  on*  reverberation?  This  report 
is  a. summary  of  the  USIj  experimental  and  theoretical  results  pertaining 
to  propagation,  reverberation  and  ambient  noise.  Published  results  of 
other  investigators  are  also  included  for  comparison. 


In  the  Arctic  Ocean,  the  speed  of  sound  is  usually  an  increasing 
function  oi  depth.  The  propagation  of  a  transient  signal  can  be  described 
by  rays* which  ^re  refracted  upward  and  suffer  repeated  reflection  at  the 
surface  (RSR) .  At  long  range  the  deepest  rays  arrive  fir st.  This  earli¬ 
est  arrival  may  be  limited  by  the  bottom,  or  by  a  "critical  angle"  of- re¬ 
flection  at  the  ice  surface.  As- time-progresses,  rays  arrive  with  de¬ 
creasing  vertex  depth  and  with  increasing  rate  until  the  train  terminates 
with  ray 8  traveling  near,  the  surface. 


If-the  first  rays  are  bottom  limited;  the  signal  onset  is  governed  by 
the  bottom,  grazing,  ray .,  For  .'greater  grazing  angles,  the  rays  are  reflec¬ 
ted  at  the  bottom  and  the  groupvelocity/is  less,  Thus  bottom  reflection 
modes  appear  as  a  sequence  of  discrete:. arrivals,  withrthe  spacing  governed 
by  grazing  angle  and  water  depth.  These  over!;  >  the  refracted-mpdes,  .and 
often  persist-for  a  longer. -time.  When  the  surface  grazing  angle  is  greater 
than  the  " critical  angle"  however;  the  bottom  modes  are  rapidly  attenuated. 
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At  long  range,  the  normal  mode  theory  also  provides  a>  simple  des- 
cription  of  the  refracted  wave  dispersion  (Sec.  2.  5) .  Because  of  the 
losses  produced  by  repeated  scattering  at  the  surface,  ordinarily  no  more 
than  two  or  three  modal1  harmonics  are  observed. 

At  short  range,,  however,  the  normal  mode  calculations  are  more 
complex  because  the  arrivals  are  less  uniformly  dispersed.  Since  the 
scattering  loss  is  also  much  less,  a  greater  number  of  harmonics  and  a 
wider  range  of  frequencies  are  received.  Convergence  phenomena  may 
also  occur  which  produce  concentrations  of  energy  periodic  in,  range. 

Thus  we  may  have  one  or  more  group  velocity  minima  in  normal  mode 
theory  corresponding  in  ray  theory  to  inflections  at' particular -surface 
grazing  angles. 

Imthe  simplified  theoretical  method  for  estimating  transmission 
loss,  we  consider  only  the  effects  of  refraction  arid  scattering  (Sec.  2. '6). 
Since  the  predominant  energy  at  moderately  iong  ranges  arrives  near  ter¬ 
mination  of  the  signal,,  spreading  loss  is  approximated  for  ,the.  pr eVailing 
half-channel  conditions,  while  the  scatteririg  loss'  is  taken  as  the  loss  per 
reflection  (Marsh)  tirnes  the  number  of  reflections.  For  the  scattering 
model,  we  use  the  wind-driven  sea  surface  model,  with  the  R.  M.  S.  rough¬ 
ness  adjusted  to  fit  the  experimental  data. 


iv  3  Reverberation 


The  validity  cf  the  scattering  model  is  confirmed  by  the  roughness-? 
wave  number  spectrum  calculated  from  measurements  of  backscatter  (re¬ 
verberation)  froth,  an  explosive  source  (Sec.  3.  2):  The  scattering  otrength 
depends  on  wavelength  and  angle  of  incidence,  so  that  the  roughness -wave¬ 
length  spectrum  levels  are  calculated,  from  the  reverberation  frequeiicy 
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spectrum.  The  procedure  is  to  determine  the  scattering  strength  as  a 
function  of  time  and  frequency.  The  surface  grazing  angle  (for  single 
incidence)  is  then  calculated,  considering  shot  depth  and  refraction  ef¬ 
fects.  The  incident  sound  wavelehg.h  and  grazing  angle  then  determine 
the  corresponding  wavelength  of  the  surface  roughness  spectrum.  The 
roughness  amplitude  is  finally  calculated^fr oro  the  scattering  strength1 
for  the  appropriate  wavelength. 

1.  4  Noise 

A  summary  of  amhient-noise,  measurements  is,  presented  in. Sec.  4.  2. 
Iix  undisturbed  periods,  the  noise  background  in  the  Arctic  Ocean  is  often 
below  sea  state  zero  equivalent.  Under  these  conditions,  measurements 
with  the  DT99  were  system-noise  limited.  For  this  reason  a  special  low 
noise  system  was  developed  for.  ambient  observation,  which  permitted 
measurement  to  some  40  dB  below  the.Knudsen  sea  state  zero  spectrum. 
For  slightly  disturbed  conditions,  the.  spectra  are  similar.  to  Knudse^s. 
During  prolonged  quiet  periods,  however,,  the  high  frequencies  drop  off 
markedly,  leaving  a  noticeable  peak  in  the  "pass  band"  of  the  sound  '.chan¬ 
nel.  This  is- interpreted  as  noise  propagated  from  distant  sources.  .Com¬ 
parison  is  also-made  between  acoustic  and  seismic  ambient  -noise'  measure¬ 
ments. 

Other  interesting  features  of  ambient  noise  unique. to  the  Arctic  in¬ 
clude  thermal  "popping!1  and  explosive  cracking  after  a  rapid  drop  in  i&r 
temperature.  Blowing  snow  noise  during  winter  storms  may  also  reach 
very  high  levels  (sea  state  3  or  more).  Generally,  however,:  the  sumriier 
calms  are  characterized  by  very  quiet  conditions.  .At4hese  times,  sounds 
of,a  great  variety  of  marine  life  are  oCcasionally  hcard.  The  sources  of 
these  whistles,  .glides,  chirps,  grunts,  etc.,,  have  not  been  identified,  but 
probably  include  cetaceans,,  seals*  walrus,  etc. 


2..  Arctic  propagation 

2.  1:  Experimental  Procedure 

In  the  propagation,  experiments  between  the  drift  stations,  shooting 
schedules  were  prearranged  or  coordinated  directly,  by- radio,  communi¬ 
cation.  Electrical  detonation  was’ used,  and  both  ifiring  and  arrival  times 
were, measured  by  chronometer  to  an  accuracy  of  -  ,0.  5  sec.  Since  the 
stations  were  continuously  drifting  in  an  unpredictable  manner,  it, was 
impossible  to  determine  the  range  between  stations  to  equal. accuracy. 

Sun  positions  were  takerinby  theodolite,  whenever  sky  conditions  permitted, 
but; the  average  error  was  no  less  than  a  few  kilometers.  Daily  drift's  on 
the  order’ of  a  magnitude  greater  complicated  the  navigation  problem. 

In  the  aircraft  runs’,  practice  depth  charges^ PDC)  were  dropped  at 
intervals  along  a  predetermined  track  whenever  suitable  openings  in  the 
ice,  cover  could  be  sighted.  The  PDC's  were  set  to  explode  at  a  depth  of 
15  or  60  meters*  and  the  detohation»tirne  was  taken  as  the  time  of  impact 
plus  the  estimated  sinking  time.  Aircraft  drop  positions  were  determined 
by  sun  shots  and  deac.  reckoning,  and  are  considerably  less  accurate  than 
the  positions  of  the  islands. 


At  the  drift  stations,  the  shot  signals  were  detected  by  means  of 'hy¬ 
drophone  slower- ed  througb-!holes  drilled  in  the  sea  ice;  Since  the  ice 
islands  are  as  much  as;  30-40  meters  thick,  a  minimum  depth  for  the  b.y- 
drophories  of  60  meters  was  chosen  to  prevent  possible  shadowing.  Some 


The  depth  settings  .of  the  PDG  pressure  actuators  are  actually  50  and  200 
feet.  For  consistency,  all  length  measurements  will  be  expressed  in  the 
metric  system. 
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measurements  were  made  with  the  phones  as  deep  as  180  meters.  How¬ 
ever,  the  greater  depths  were  found  to  be  of  marginal  interest  since  the 
simpler  procedure  of  varying  shot  depth  showed  most  of  the  essential  fea¬ 
tures. 


In  all  the  earlier  experiments,  a  standard  PQM  (DT  99  )  hydrophone 
was  used.  The  disadvantages  of  this  hydrophone  were  that  the  multi¬ 
conductor  cable  was  somewhat  unwieldy  and  the  associated  preamplifier 
was  subject  to  saturation  .by  "strong  signals.  The  overall  frequency  re¬ 
sponse  was  10-10, <000  Hz.  For  later  teats,  a  new  system  was  developed 
using  a  bilaminar -disc  hydrophone.  Because  of  their  high<  capacitance 
(O.lpF) ,  these  units  could,  efficiently 'drive  several  hundred  feet  of  light¬ 
weight  shielded  cable  (Hrectly  without  a  preamplifier.  The  system  noise 
was  much  lower  than  it  was  using  the  DT99;  however,  high  frequency  re¬ 
sponse  was  sacrificed  since  resonance  occurred  at  2  kHz.  Another  disad¬ 
vantage  was  that  the  hydrophone  xould  not  withstand,  more  than  100  meters 
submergence  without  mechanical-failure. 

The  signals  from  the  hydrophone  were  recorded  by  means  of  a  graphic 
recorder,  with.fr equency  response  to  100  Hz.  Magnetic  tape  recordings 
(  FM  and/or  direct)  were  also  made  with  full  system,  frequency  response  for 
later  analysis.  Calibration;  of  the  system  was  accomplished  by  injecting 
sinusoidal  signals  of  known  voltage  at  selected  frequencies.  A  single  master 
gain  control  provided  adjustment  of  the  system  gain  for  calibration,  signal 
and  ambient:  noise  recording; 

During  the <earlier  tests,  it  was  found  that  the  dxplesions  often  de¬ 
livered  varying  yields,  as  result  of  incomplete  detonation.  To  circumvent 
this  difficulty,  a  procedure  was  established  to  record ,th(5  shot  at  the  trans¬ 
mitting  site  whenpoasibVe,  The  effective  yield  was  determined1  from  the 
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bubble  pulse  period  and' depth.  The  record  in  Figure  2  is  a  pressure  vs 
time  signature  for  a  typical  shot. at  1  km  showing  the  shockwave  and  bubble 
pulses  (positive  pressure)  ,  and  surface  reflections  (negative  pressure)  . 

If  the  system. were  pverloaded  or  cluttered  with  interfering, arrivals,  the 
bubble  pulse  period  could  be  obtained  from  the  periodic  envelope  modular 
tion  Of  the  spectrum.  The  two  spectra  shown  in  Figures*  3a,  b,  which  are 
shot  signatures  at  1  km' and  50  km  range  respectively,  show  this  clearly. 

2  '2  Data  Analysis 

Analysis  of  the  data  is  bf  three  general  types.  First,  graphic  re^ 
cordirigs  made  on  the.  spot  furnished-direct  pressure  amplitude  vs  time 
traces.  These  show  qualitatively  the  effects,  of  shot  size  and  depth,  as 
well  as  the  qua  si- sinusoidal  and  interference  features  of  the  propagation. 

The  times^of  onsetcand  termination  also  provide,  measures  of  time -disper¬ 
sion  and'  travel  time. 

Secondly,  the  magnetic  tape  recordings  were  analyzed  by  means  of 
a  Kay  Mis  silyzer  Spectrum- Analyzer  to  provide  the  frequency  spectrum 
v.s  time.  These,  graphs  are  useful  in  analyzing  the  frequency -time  depen¬ 
dence  of  the  normal  modes  and  harmonic  relations  betweenthem  as,  well 
as  the  interference  effect  of  the  bubble  pulse; 

Finally,  the  magnetic  tapes  were  used  to  make  continuous  loop  re¬ 
cordings  of  the  shot  signals;  which  were  scanned  by  means  of  a  variable 
frequency,  narrow  band  filter.  A  typical  result  of  this  technique  is  shown 
in,  Figure  4.  Together  with  the  calibration  tapes,  this  permits  calculation 
of  the'  received  energy  flux  spectrum  as  a  function  of  range,  shot  size, 
depth,  etc.  Knowing  the  source  energy  flux  spectrum  at  1  m,  range-.(Weston)5, 
the  transmis  sionlos  si's  calculated.  Since  the  received  signal  varies  in 
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time  and  frequency  in  a  complex  fashion  because  of  source  and  modal  in¬ 
terference  effects,  the  peak  envelope  method  of  analysis  was  used' in  much 
of  this  work.  A  smooth  curve  was  fitted  tO'the  "average"  spectral  peaks 
so  as  to  eliminate  mode  and  bubble  pulse  modulation.  This  is'the  solid 
curve  seen  in  Figure  4.  In  another  technique,  logit  filter sL  were  used/ and- 
the  total  (integrated)  energy  flux  spectrum  levels  calculated.  The  two 
methods  were  found  to  be  in  quite  reasonable  agreement. 


2.  3  Experimental  Results  (General); 

The  pertinent  data  for  the  -analyzed  shots  are  contained  in  Table  I. 
Each  record  is  identified  both  serially  and  with  the  notation  of  the- origin 
nal  experiment  log. 

Figure  5  shows1  typical  graphic  recordings  made  at  ARLISII  from 
shots  of  several  sizes.detonated  af/T-3.,  All  were  "at;  120  meters , depth 
except  the  375  lb.  shot,  which  was  fired  at  100  meters.  The  travel  time, 
corresponding  to  termination  of  the  water  wave  is,  563  -sec.  For, -a  sound 
speed  of  1/438  m./sec,  the  calculated  range  is:  810  km,  agreeing  with  the 
navigational  value  within  a  few  km.  It  should  be  noted  that  the  receiving 
system  gain  was  adjusted  to  compensate  for  source  level  changes.  The 
timecdispersion  varies,  between  6.  6  and  7.  1  seconds,,  increasing. with  shot 
size.  For  the  smaller  shots,  the  complexity  of  the  beginning  of  the  water 
wave  is  particularly  noticeable.  This  is  due  to  contributions  from  several 
modes  as  well  as  from  interference  of  the  bubble  pulse.  For  the  larger 
shots,  the  spectral  energy  is  concentrated  at  lowex*  frequencies,  the  rela¬ 
tive  excitation  of  higher  frequencies  is  less  and, the  bubble  puis e  frequency 
is  lowers  Source  interference  effects  persisting  throughout  the  entire  train 
are  clearly  evident. 


Figure  6  shows  Missilyzer  recordings  of. two  typical  shots.  Note 
that  the  lowest:C  zero)  mode  begins  neaTr  10  Kz  and  increases  with  time  to 
about  70  Hz  at  termination.  The  first  modal  harmonic  is  visible  through¬ 
out  the  trace,  while  the  second  and  higher  modes  vanish  soonafter  the  be¬ 
ginning  of  the  wave.  Bubble'  pulse  interference  effects  are  evident  from 
the  amplitude  modulation  of  the  mode  lilies  at  multiples  of  the  bubble  pulse 
frequency.  Norma!  mode  theory  is  discussed  in  Section  2.  5. 

The  received  signals  characteristically  have  an  abrupt  beginning  and 
end  with. a  dur atiori  pr opor  tional  to  travel  time.  Figure  7  contains  a  plot 
of.  time  duration  (dispersion) Vs.  travel  time  for  some  fifty  shots,  detonated 
under' various  conditions  and  times  during  the  experimental  program.,  The 
data  show  a  mean  time  dispersion  ratio  of  slightly  rpore  thah  1%.  The  fac¬ 
tors  governing  onset  and  termination  of  the  water  wave  are  discussed  else¬ 
where; 


Figure  8  give 8  a  typical  sound> velocity  "profile  for  simnmer  conditions5 
The  two  major  zones  of  interest, ar.e  the  deep  water  isothermal  gradient  be? 
ginning  near  300  meters  and  the  overlying  layer  of  some  3  or  4  times  greater 
average  gradient.  In  the  upper  zone,  there  is  some  detailed  structure  which 
depends  on  season. 

Table  II  contains  ray  calculations  for  tie  velocity  profile  of  Figure  8 
(see, Section  2.  5).  It  is  notable  that  the  time  dispersion  corresponds  to 
ray  depths  less  than  1500  meters  and  surface  grazing  angles  less' than  13°. 
Since  the  minimum  water  depths  for  some  of  the  experimental  paths  were 
in  excess  of  2000  meters  (Figures  9  —  1 1 )  ,  it  is  clear  that  the  onset  of  the- 
wave  train  is  not  always  governed  by  the  bottom,  If  thi  s  is  true,  the  on¬ 
set  must  be  limited  by  the  surface-  reflection  coefficient  of  the  ice  cover, 
which  apparently  diminishes  rapidly  for  grazing  angles  greater  than  13° 
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Figures  12a,  b  show  Missilyz'er  traces  tor  the  propagation  path.pro- 
files  of  Figures  11a,  c.  In  lie,,  anunderwater  obstacle  {the  Chukchi  Rise) 
shadows  all  rays  deeper  than/'20&  meters,  leaving  only  one  second  of  the 
dispersion  pattern  remaining.  The  deep  water  duration's  in  excess  of  ten 
seconds. 

Under  some  circumstances,  the  bottom  reflection  mode  is  received 
along  with  the  RSR  mode;  Lythis  caiae  the  onset  corresponds  to  a  group 
velocity  maximum.  This  effect  .is  seen  in  Figure  13a  which  is  a  recording; 
of’ a  shot  dropped  from  an/air craft:  and  received  572  seconds  later.  The 
first  hundred  kilometers  of  the  propagation  path  were  in  shoaling  water  (10a)' 
and  gave  rise  to  .more  than  twenty  identifiable  bottom  arrivals  which;  obscure 
the  RSR  mode.  Figure  13b  ir  another  shot  at  almost  the  same  range  but  re¬ 
ceived  some  20  minutes  later.  The  bottom  topography  { 10b)  has  changed, 
so  that  an  obstacle  intervenes  and  cutsoutthe  bottom  mode  completely: 


2.  4  Transmission  Loss  { Experimental;) 

Figure  14  shows  a  typical  energy  flux  vs., frequency  spectrum  (L>  ) 

•  ’  '  ’  ■■  r 

measured  as  described  in  Section  :2.  2.  Transmission  loss  vs.  range  data 

for  all  shots  tcomputedr from  the  relation, N  =  L  -  h  Yare  shown  in 

y/  s  r 

Figures  15  through  27.  The  source  level  s.oectr.um  at  1  meter,  L  ,  was 

"  r  ' 

taken  from  Weston  for  the  appropriate  yield,  correct ed/( wh e re,  po a s ibl e) 
by  the  bubble-pulse-peripd-ys-depth  method  mentioned  earlier:  It  should 
be  noted  that  Weston's  shot  spectra  account  for  bcbfclc  pulse  interference 
only  at  the  fundamental  frequency,  while  the  true  spectra  (shows  maxima 
and  minima  over  many  harmonics.  Thus,  the  Weston  spectrum. represents, 
a  quasi-envelope  of  the  actual  one  (Figure  28)  .  The  harrow  band  analysis 
method  used  here  clearly  shows  these  ana  other  interferences  and  an  en¬ 
velope  was  drawn  similar  to  Weston's  t&average  them  out. 


The  transmission  loss  data  show  considerable  scatter  which  is 
taken.to  represent  the  combined  effects  of  variations  in  the  transmission 
path  and  ice  cover,  asweli  as  experimental  error.  Some  of  the  short 
range  levels  are  considerably  higher  than  predicted  for  the  half  channel 
and  are  interpreted  as  possible  results  0/  convergence.  This  effect  is 
'illustrated  in  Figure  r29  a,  b  where  arrivals  from  #8  caps  at' two  different 
ranges  are  shown.  The  first  at  1  km  shows  only  a  discrete  arrival  fol- 
lowed  by  a  bottom  reflection.  The  second  at  5  km  also*  apparently  shows 
a  delayed  convergence  arrival  for  which  the  level  is  about  20  dB  higher 
than  the, direct  arrival.  Buck6  has  also  observed  high  signal  levels  at 
short  range  which  may  be  ascribed  to  convergence. 

Milne! s7  transmission  losses  out'to  90 -km  have  been  included  in 
the  figures  for  comparison.  These  measurements  were  made  in  the  spring 
and  .may  represent  differences' in  the  ice  condition*  or  a  change  in.the  re¬ 
fraction. 

The  theoretical  transmission  loss  curves  are  discussed  in'Section 
2.  6.  The  curve  'S'  is, the  surface  scattering model  representing  the  "aver? 
age"  loss  for.  a-fong  transmission  path.  The  curve  'A1-  is  absorption  loss 
only '.and  represents  a  lower  limit  for  short  range  where  no  surface  reflec¬ 
tions  are  involved; 

Various  factors  may  be  .responsible  for  the  transmission  loss  "anoma¬ 
lies".  For  example,  in  deep  water  the  early  arrivals  may  dominate  and  the 
high  frequency  loss  may  be  low  because<the  number  of. reflections  is  less. 
Thus  the  transmission  loss  depends  on  the  bottom  topography  as  well  as  the 
condition o.f  the  surface  and  the  sound  velocity  profile. 

However,  the  scatter  of  the -transmission  loss  data  is  evidently  not  due 
to  short  term  changes  ip  the  medium,  which,  in  fact,  is  remarkably  stable. 
Figure  30  shows  a  series  of  f  ib.  shots  at  80  meters  depth  received  over 
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an  800  kilometer  path.  The  two  hydrophones  (a,  b)  were  1  km. apart  and 
the  time  between  shots  was  2  hours.  The  similarities  indicate  both  space 
and  time  stability  of  the  acoustic  path. 
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,2.  5  Normal  Mode  and  Ray  Theories  8 


To  illustrate  the  .connection  between  normal  mode  and  ray  theories,, 
we  consider  firct  a  single  positive  .gradient  defined  by  c  =  c  0  (T  +  gz^}.,. 
where  c_  is  the  depth  dependent  sound  velocity,  g  is  the  velocity  gradient 
(ni  ),  c  is  the  sound  velocity;  (ms”  )  at  the  pressure  release  surface 
and  z  is  the  depth  (m) .  The  resulting  separated  equation  for  the  depth 
dependent  portion  of  the  velocity  potential  is  . then  approximately , -given  by: 


+  ^<i>  =  0 


where 

X  = 

rv>  jv 

\  c„  2 

/  tcT* 

(1) 

(2) 


Here,  w  is  the  angular  frequency  (  s"  1 ) ,  and  K  is  the  wavenumber  (m”  1  ) 
The  eigenfunctions  of.  Equation  (  p)  are  modified  Hahkel  functions  of  order 
V3  •  hj  (  (, )  and  h2  (  £,)  represent  "downward"  and  "upward"  going 
waves.  For  our  purposes  we  consider  o.nly  ht  (  £  ) .  The  eigenvalues  are 
obtained  fromthe  zeroes  which  have  the  argument  7r/3  \ Figure  31). 


Thus,  for  the  pressure  release  surface,  the  first  mode  4>  =  0, 
z  =  0,  n  =  1  is  identified  .as  point  (  a) ,  while  the  curve  to  the  right  repre*» 
sents  the  depth/ frequency  dependence  which  increases  to  its  maximum  and 
then  decays  exponentially.  Similarly,  the  curve  to  the  right  of  point  (b) 
represents  the  second  mode  n  =  2. 
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Good  approximation  of  the  zeroes  of  hj  are  given  by 


3  7T 


2  (n  -  |  )  n  =  1,2,3 


(3) 


th 


where  x  =  E,  ( z  =  0).  From  this,  the  phase  velocity  of  the  n  mode 
becomes  z 


k 


=»*£  -  <*.(*♦+  I <» -i »]  ) 


(4) 


and  the  .group  velocity,  u^,  is  obtained  from  the  relation 

d'c 


n 

u  =  c  +  u  — — - 
n  n  dw 


(5) 


to 


Substituting  f  =  the  final  result  is.: 


u  =  c, 
n 


•  (1+i-  (n-i)]  ) 


(6) 


The  frequency-time  behavior  of  a  wave  received  at  a  distant  sur¬ 
face  point.from  a  surface  impulsive  source  is  simply,  the  frequency -mode 
number  components  delayed  according  to -the  appropriate  value  of  the 
group  velocity. 


To  compare  ray  theory,  the  familiar- equations  for  the -positive 
gradient  g  is:, 


To 


2  ,  sin,  0O 
g  cos  :0O 


(7) 


where  r0  is  the  horizontal  distance  between  surface  reflections  and  0 
is  the  grazing  angle.  Also, 

0o 

t  2  f  de  ~  Ju 

0  ~  c0  g  *  CO®  6  =  Cog 

0 
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where  t0  is  the  travel  time  for  the  distance  r0.  From  this,  the  group 
velocity  u  is  obtained: 

u  =  c°  (  1  +  ^  *  *  *  *  )  (9 ) 

To  determine  the  frequency  f,  the  reciprocal  of  the  time  T  between 
the  Nth  and  {N  +  l).th  arrivals  at  a  distance  R  is  taken.  That  is: 


f 


dN 

dT 


(10) 


Since  the  range  R  =  Nr0  and  the  travel  time  X'  =  Nt0  ,  then 


- 1 


f  &  to 


du 

dr0 


(.11) 


From  the  above, 


u  =  c0 


( >  *  f  ) 


and  03  = 


3  cp  g 

Z£ 


(12) 


Combining  the  previous  results  yields 


73 


*  a  c«  (1+  t  ’) 


(13) 


Comparing  Equations  (6)  and  (13)  ,  we  find  the  only  difference  be¬ 
tween  the  two  lies  ,in  the  mode  number  term.  In  fact,  the  ray  frequency 
is  3/4  the  lowest  mode  frequency. 


For  the  bilinear  case, 


c  =  C0(l  +  g07.) 


cz  =  Cj  (1  +  g.  2,/ 


0  <  z  <  zi 
Zi<  z  <  ©• 


(14) 
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The  approximate  equation  analogous  to  Equation  3  is: 


/  V 

/?-  +  a(x  -  1)  /z  =  b 


where 


=  (go  Ci2  / gi  c02  -  1) 


I  1 

(n  -T  '  ,  TT 
(  2  So  zi  )  2 


The  resulting  group  velocity  is: 


u  S.  Cq  [  1  +  g0  zx  (x  +  "  )] 


=  C0  T}-+,gO'zl 


bA  +  q(  x-l)  S  1 

^  +  «  U-1>  *  j 


In  the  above,  x  -  1  is  replaced  by  zero  when  x  <  1,.  and  in  this  case 

'•  (iti  !¥  ‘t-i'f)  1 


whfch  is  the  same  as  Equation  (  6)  . 


A  graphical  solution  of  Equation.  15  is.  readily  obtained  with  the  aid 
of  the  nomographs  (Figures  32, ,  33).  To  obtain  x,  use  the  auxiliary 
equation  b  =  cr/f,  locate  f  on^the  abscissa  of  Figure  32,  prqceed-  verti¬ 
cally  to  the  appropriate  b  curve  (marked  with  the  value  a)  ,  proceed  hori¬ 
zontally  to  the  appropriate  a  curve,  and  thence  downward  to  read  x  on 
the  abscissa;,  uiyen  x,  locate  the  abscissa  in  Figure  33,  obtain  a (x;-l)  2 
as  the  corresponding  prdinate,  locate  b/3  as  the  ordinate  corresponding, 
to  the,  chooen  f  (abscissa)1;  obtain  x  2  +  qlx  -  l)  2  as  ordinate  corres¬ 
ponding  to  abscissa  x,  and  complete  the  calculation.  Figure  34  shows 
the  resulting,  frequency  dispersion  vs.  group  velocity  for  this  bilinear 
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gradient  model  with  appropriate  values  calculated  from  Equation  16. 

The  corresponding  ray  equations  are: 

r°  *  g0"cos'e"  .[(*&! )  sin-Oi  +  sin  0O  -  sin  0tJ  (18) 

t0  =  — —  f  (  ln(l  +  sin  0!  )  -  ln(l  -  sin  0,  ) 

go  co  w.  gi 

+  In  {1  +  sin  0O )  -  In  (T  -  sin  0O)  )]<») 

The  calculations  for  the  multilinear'  gradient  (Table  -IT)  were  ob¬ 
tained  by  a  ray  computer  program  for  ten  gradients  fitted. to  the  depth 
profile  of  Figure  8;  The  frequencies  were  calculated  from  Equation'(ll) . 


Transmission  Los3  Theor\ 


Earlier,  the  bottom  grazing  ray  for  the  bilinear  gradient  was  iden¬ 
tified  with  a  group5 velocity  maximum.  This  determines  the  time  of  onset 
of  the  wave  train,  provided^the  surface  grazing/ angle  is.  less  than  13°  .  For 
the=  bilinear  gradient  (large  -  small  )?model,  this  is  the  only  group  velocity 
extremum,  except  for  the  trivial  minimum  at  termination'  of  the,  wave  where 
the' ray  vertex  depth  goes  to  zero  and  the  frequency  incr eases/ without  limit. 
In  the  multilinear  gradient  case,  the  situation  is  quite  different  since  there 
is  a.  group  velocity  minimum  for  a.  surface  grazing  angle  of  7.  7°.  Since 
the  vertex  depth  for  this  ray  is  only  about  240  meters,  no  significant  con¬ 
vergence  is  expected  if  source  or  receiver  are  deeper  than  this.  ILthey  are 
less  deep  the  main  effect  at  long  range  wilibe  a  sudden  and  apparently  pre¬ 
mature  termination  of  the  wave  train.  At  short  ranges,  however,  conver¬ 
gence  focusing  is  important.  Thus,  as'the  range  increases,  we  expect  con¬ 
centration  of  energy  at  multiples  of  the  skip  distance.  After  a  few  surface 
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bounces,  .zones  should' tend  to  spread  and  overlap  because  of  scattering, 
so  that  the  decay  eventually  becomes  a  monotonic  function  of  range.  At 
long  range,  therefore,  near  the  termination  of  the  wave  train,  the  total 
transmission  loss  may  be  estimated  from  the  combined  effects  of  refrac¬ 
tion  and  scattering  as  nearly: 

Nw  =  10  log  r0  +  10  log  R  +  Nag  (20). 

where  r0  is  the  ray  skip  distance, R  is  the  range  in  meters,  N  is  the 
number  of  surface  incidences  and  ag  is,  the  loss  per  reflection.  As  a 
scattering  model;  the  wind  generated  ocean  spectrum  is  taken  to  ap¬ 
proximate  the  ice  roughness.  According  to.  Marsh? 

„4  3/  8/ 

a  =-  10  log  (1,-  2.  57  X  10  f  /2h/5  sin  e  )  (21) 

s 

where  f  is  the  frequency  (Hz):and  h  is  the  R...M.  S.  roughness  in  meters. 
For  a  small  loss  per  Bounce,  (20 ^becomes  approximately: 

3/8/ 

Nw  S  20  log  r0  +  10  log  N  +  N  (1.  10  x.KT  3  f  '  z  h  !s  sin  0  )  (22) 

where  R  =  Nr0 

To  give  specific  results,  take  0  =  7.  7° ,  r0  =  6021.meters, 
h  =  2.  4:  meters,  and  N  =  134,  Then 

7  3  / 

N  »’W'+Rx.l0"  f  /z  (23) 

w 

For  these  conditions,  which  correspond  to  shot  31B  at  810  kilometers, 
the  following  is  calculated: 
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Figure  35  shows  a  plot  of  Equation  ( 23 )(  solid  line  ), :  together  with  the  ex¬ 
perimental  transmission  loss  values  (circles)  computed  by  N  =L  -  L 

w  s,  r 

(Section  2-4).  In  the  frequency  range  shown,  the' fit  is  quite  good;  however, 
above  100  Hz,  the  theoretical  scattering- loss  apparently  becomes  excessive. 

Figures  15  through  27  show  the  compiled  transmission  loss  data.  In 
Figures  15-24,  the  curve  S  is  Equation  ( 22) ,  The  straight  lines  are  the 
first  two  terms  of  Equation  20,  giving  the  spreading  loss  only.  The  usual 
absorption  loss  is  the  curve  A  in  Figures  24-27.( which  includes  scattering 
loss) .  For  the  water  temperature  of  0°  C,  a-bsorption  loss  is  taken, as 
012fzdB/km,  where  f  is  the  frequency  in  kHz.  The  curve  A  is  included  to 
show  the  lower  limit  of  transmission  loss  in  the  absence  of  scattering  (high 
frequency  -  short  range).  Both  the  5  and  A  curves  are  shown, in  Figure  24 
to  illustrate  the  transition  region. 

Equation  22  appears  to  give  a  fairly  pessimistic^approximatioriand  is 
probably  applicable  when  the  water  is  le.ss  than  500  meters  deep  over  a  por¬ 
tion  of  the  transmission  path.  For  deeper  water,  the  high  frequencies  tend 
to  be  dominated  by  the  deeper  rays  (earlier  arrivals)  which  perhaps  accounts 
for  the  low  loss  anomalies.  It  is  doubtfvl,  however,  that  more  detailed, cbn- 
sideratibn  of  ray  paths  and  bottom  toppgraphy  would  significantly  improve 
the  picture  at  this  point,  since  the  grazing  angle  dependence  of  the  surface 
is  not  known. 


3.  REVERBERATION 
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3.  1  Experimental  Results 

The  results  reported  here  concern  the  local  scattering  or  reverbera¬ 
tion  produced' by  the  rough  under  side  of  the  ice  in.the  vicinity  of  explosive 
sources.  The  explosives  were  monitored  and  recorded  through  hydro¬ 
phones  some  hundreds  of ’meters  away.  Table  HI  gives  the  pertinent,  data 
identifying  the  shots  and  recordings.  The  magnetic  tape  records  were  sub¬ 
sequently  played  back  through  a  Bruel-Kjaer  analyzer,  using  l/3-octave 
(10%  bandwidth)  filters -centered  at  frequencies  between  40  and  8000  Hz. 

The  filtered  waveforms  of  shot  8  and  a  composite  of  shots  10-12  are 
shown  in  Figures  36  and  37,  respectively.  It  should  be  noted  that  the  time 
scales  of  the  two  differ  by,  a  factor  of  10.  In  the  figures,  time  is  marked 
fronvthe.  shot  instant.  Figure  37  shows  an  early  arrival,  particularly  at 
the  lower,  frequencies,  which  is  probably  an  ice  wave;  Shortly  after  that 
comes  the  direct  and  surface  scattered  wave,  which-decays  for  a  second 
or  so  before  the  first  bottom  arrival  is  seen.  Additional  orders  of  bottom 
reflections. may  be  seen  a:  bit  later  in  the  figure.  In  Figure  36,  these  arri¬ 
vals  are  compressed  within  the  first  few  seconds  of  the  record,  and  the 
main  curve,  lasting  up  to  nearly  two  minutes,  consists  of  reverberations 
returning  from  the  ice  and  the  bottom. 

Quantitative  results  are  conveniently  presented  in  terms  of  the  scat¬ 
tering  strength,  Ng. 

Ng  =  -  Lg  +  2N^.^  -  10  log  A  (24  ) 
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In  Equation  (24) ,  is  the  level  of  reverberation  at  time  t  (  dB/^ibar  )>, 
bc  is  the  level  of  the  source  (dB//pbar.)  ,  2N  is  the  propagation  loss 
from  source  to  reverberating  area  plus  loss  from  reverberating. area  to 
receiver,  and  .A  is  the  effective  reverberating  area  at  time  t.  The  .re¬ 
verberating  area  is  the  elliptical  annulus:  on  the  surface -of  width  deter¬ 
mined  by  the  duration  of  the  incident  wave  and  generally  located  by  the 
travel  time  from  source  to  surface  to  receiver.  Table  IV  presents  se¬ 
lected  values  of  Ng  tabulated  from  the  shots  of  Table  HI,  Marsh  has 
shown  how  the  scattering  strength, may  be  expected  to  depend  upon  the 
properties  of  the  scattering  surface,  and  upon  the  geometry  of  the  situa¬ 
tion.  With  the  bistatic  setup  prevailing  in  the  present  experiments,:  it  is 
quite  difficult  to  interpret  directly  the  scattering  strength  for  comparison^ 
with  other  geometry  or  other  surfaces.  For:this  reason,  the  ice  rough¬ 
ness  spectra  is  derived  from  the  observed  scattering  strengths,  and  then 
the  scattering  calculated  for  other  geometries.  This' is  done  below. 


For  L  the  value s/o£  Weston  are  used  of  the,  spectrum  level  for  the 

\V> 

appropriate  yield,  to  which  20  log  (bandwidth)  is  added.  The  yields  were 
determined  experimentally  from  missilyzer  recordings,  in  which  the 
bubble  periods  were  evident.  Theue  yields  in  several  cases  differ  marked¬ 
ly  from  the  design  yield  based  on  weight  of  explosive.  2N^^  is  taken  , to 
be  simply  40  log  ct,  c  being  the  mean  velocity  in  meters  per  second  Fi¬ 
nally,  the  reverberating  area  is  approximately  xc2t  T  / 2,  X  beirig  the  ef¬ 
fective  pulse  length,  twice  the  reciprocal  bandwidth.  These  approximations 
are  adequate  to  rectify  the  experimental  values  of ;Ng,  producing  a  quantity 
which  is  a  function  only  of  fr<  juency  at  a  given  time,  asd whose  variation 
with  time  can  be  interpreted  in  terms  of  the  prevailing  propagation  condi¬ 
tions. 


20 


3.  2  Reverberation  Theory  and  Interpretation 

Marsh9  gives  the  scattered  intensity  I  in  terms  of  the  intensity  I  of 
waves  incident  upon  a  rough  surface  and  the  wave  number  spectrum  of 
surface  roughness- A2  (k)  in  the  form 

1  =  JfirdXd\L 

T  =  y2k3A2  (ka) 

The  direction  cosines  pf  the  incident  wave  are.oc,  (3,  y,  and  those  of  the 
scattered  wave  are  X,  p,  v  .  So,  fpr  instance,  the  sine  of  the  grazing 
angle  of  incidence  is  y  =  sin  0  .  The  quantity  q  is:  [(*-«>2+(p-p)2]  2 

The  reverberation  returning  at  times' greater  than  about  10  seconds 
comes  from  an  approximately  circular  annulus..  The  separation  of  source 
and  receiver' may,  be' ignored.  Under  these  conditions,  the  integral-may 
be  carried  out  explicitly.  In  so  doing,,  effects  of  refraction  must  be  con¬ 
sidered;  For  the  transmission  times  involved;  the  law  of  propagation 
may  be  taken  as  approximately  inverse  square.  However,  refraction  has 

x  i.  *  , 

a  strong  effect  upon  the  grazing,  angles. 

To  complete  the  calculation,  it  is  necessary  to  relate  the  surface 
grazing  angle  to  time.  For  this'  purpose,  the  bilinear  gradient  profile  ap¬ 
proximation  is  employed.  Figure  38  displays  the  surface  grazing  angle  as 
a  function  of  time  after  the  first  arrival.  Following  Mar sh  for  the  reduc¬ 
tion  of  the  integral,  we  finally  obtain,  approximating  cos  9  -by  1, 

10  log  A2  (2k)  -  N  -  30  log  k  -  40  log  sin  9+3  (27) 

8 


(25) 

(26) 
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Values  of  A2  (k)  reduced  using  this  equation  are  shown  in  Figure  39. 
Those  values  for  k  =  0.  35  and  0.  60  are  obtained  from  shots  10  -.  12, 
during  the  time  preceding  the  first  bottom  arrival.  For  these  values, 
the  bistatic  geometry  is  important,  and  the.  integration-corresponding 
to  Equation  (25)  cannot  be  carried  out  in  closed  form.  However,  the 
variation  of  N  with  frequency  would  be  nearly  independent  of  the  bi- 
static  geometry.  Accordingly,  these  values  have  been  adjusted  to  give 
values  of  A2  (k)  equal  to  those  of  the  other  shots  at  k  a  5.  9.  The  other 
values  for  these  small  values  of  k  are  not  plotted*  because  they  are  be¬ 
lieved  to  be  due  primarily  to  bottom  reverberation.  However,  the  cor¬ 
responding  values  of  scattering  strength  are  included  in  Table  IV. 


It  may  be  noted  that  the  rms  roughness  er  of  the  surface  is  given 


by 


<r2  =  A2  (k)  dk 

J  o 


(28) 


By  numerical  integration,  <r  pi  3m,  slightly  larger  than  the  value  derived 
from  forvvard  scatter;  but,  both  are  quite  compatible  with  ice-thickness 
measurements  of  Lyoris10  shown  in  Table  V. 


For  comparison,  the  sea-surface  spectrum  of  Burling  is  also  shown 
in  Figure  39.  The  curves  are  roughly  parallel  above  k  =  5,  and  it  appears 
reasonable  that  the  ice  spectrum  is  asymptotic  to  k"3.  This  corresponds 
to  reverberation  which  is  independent  of  frequency,  similar  to  that  in-the 
open  sea,  but  about  40  dB  higher,  A  compilation  cf  various  geophysical 
spectra  shown  in  Figure  40  indicate  -that  the  K~  3  spectrum  is  typical  of 
naturally  occurring  rpugh  surfaces.1* 


There  is  considerable  scatter  in  the  data  reported  (-standard  devia¬ 
tion  ~  5  qB)  ,  some  of  which  might  be  eliminated  by  a  detailed  computation 
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of  ray  paths  and  losses,  taking  into  account  the  actual  shot  and  hydro¬ 
phone  depths,  in  which  the  refraction  could  play  a  part. 

An  additional  feat  ure  of  the  recordings  is  the  peak  in  levels  occur¬ 
ring  from  40-50  seconds,  particularly  evident  at  the  higher  frequencies. 
This  peaking  is  undoubtedly  a  refraction  effect,  and  has  been  systemati¬ 
cally  observed  on  all  shots  identified  in  Table  HI.  The  travel  time  of  50 
second*  corresponds  to  a  grazing  angle  of  16*  for  the  direct  ray  and  an 
angle  of  11*  for  the  singly  reflected  ray.  These  angles  are  on  either 
side  of  the  "critical  angle".  There  is  thus  additional  evidence  that  re¬ 
flection  is  quite  poor  at  grazing  angles  exceeding  some  13*. 

More  recently,  other  investigator  slz,13have  made  reverberation 
strength  measurements  on. young  sea  ice  and  summer  polar  ice  over  a 
considerably  wider  range  of  grazing  angles  and  using  improved  geo¬ 
metry  to  eliminate  bistatic  problems  and  bottom  echoes.  While  all  mea¬ 
surements  agree  quite  well  foif  grazing  angles  ~  10*  (Figure  41),  the 
angular  .dependence  of  the  new  data  is  much  less  than  predicted  by  Equa¬ 
tion  (27).  This  may  not  be  surprising  since  Brown’s  normal  reflec¬ 
tivity  is  some  15  d.B  below  perfect  at  5-10  kHz. 

Kuo14  has  considered  the  Brown  experiment  by  extending  Marsh’s 
scattering  formula  to  the  imperfect  reflector  and  fitting  his  modified  re¬ 
flectivity  function  of  grazingiiangle  to  the  data.  The  results  are  shown  in 
Figure  42  and  the  derived  roughness  spectrum  in  Figure  43.  For  large  k, 
(She  agreement -between  the  Brown  and  Mellen-Marsh  roug.  cs  spectrum 
levels  is  fairly  good  even  though  Brown’s  total  roughness  is  apparently 
much  less. 
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Other  discrepancies  are  evident  from  the  work  of  Milne15  on  the 
spring  Arctic  pack  ice  shown  in  Figure  44.  Comparing  values  at  10° 
grazing  angle,  his  scattering  strengths  are  some  20  dB  higher  and  show 
even  different  grazing  angle  di  pendence.  Milne's  normal  reflectivity 
is  almost  perfect  to  1  kHz  and  only  5  dB  down  at  5  kHz  (Figure  45)  . 


4.  AMBIENT  NOISE 


4.  1  Experimental  Procedures 

Early  attempts  to  measure  the  ambient  noise  background  with  the 
DT99  hydrophone  were  usually  frustrated  by  excessive  system  noise. 

For  this  reason  bilaminar-disc  hydrophones  (USL  Type  XU  12,69  and  XU 
1333)  were  used  in  conjunction  with  a  low  noise  vacuum  tube  preampli¬ 
fier.  Since  the  capacitance  of  the  units  was  approximately  0.  1^.F  and 
the  depth  of  submergence  only  60  meters,  the  hydrophone  was  used  to 
drive  a  suitable  length  of  transmission  line  directly  and  the  preamplifier 
was  installed  above  the  surface,  often  in  the  instrumentation  hut. 

The  noise  curves  for  typical  DT99  and  XU  systems  are  shown  in 
Figure  46'which  also  includes  the,  Knudsen  sea  state  zero  curve  for  a 
refer'ence.  The-;main  disadvantage  of  the  XU  system  was  the  low  reso¬ 
nant  frequency  which  prevented  measurements  over  as  wide  a  frequency 
rangeas  the  DT99.  However,  since  the -main  concern  here  was  for  the 
frequencies  less  than  i  kHz,  this  was  not  considered  serious. 

Recordings,  of  ambient  noise  together  with  approximate  calibra- 
tionscwere  made  at  intervals.  These  recordings  were  subsequently  ana¬ 
lyzed  by  means  of  a  WE  3A4A  spectrum  analyzer  and  the  corresponding 
ambient  noise  spectra  calculated. 

The  measurements  with  the  XU  1269  were  made  at  ARLIS  H  during 
the  period  ofvSeptember  —October  1961.  Other  measurements  were  made 
between  May  -  September  1962  at  both  T  3  and  ARLIS  EU  The  data  in  no 


way  represent  a  complete  statistical  study  of  ambient  noise  even  oyer 
those  restricted  periods  because  of  limitations  set  by  the  necessity  of 
almost  complete  shutdown  of  generators  and  other  camp  activity.  Usu¬ 
ally  a  small  generator,  suitably  shock  mounted,  was  used  to  power  tht 
equipment.  Measurements  ware  concentrated  in  periods  of  light  winds; 
High  wind  speeds  tended  to  increase  water  currents  relative  to  the  ice 
causing  hydrophone  flutter  which  made  ambient  measurements  impos¬ 
sible.  Later,  cable  fairing  was  employed  to  solve  this  problem. 


4.  2  Experimental  Results 

Unusually  quiet  conditions  prevail  during  the  summer  months  when 
temperatures  are  stable  and  winds  moderate  for  a  considerable  part  of 
the  time.  This  is  in  part  due  to  the  mushy,  broken  condition  of  the  ice 
cover/ which  makes  it  a  poof  noise- generator.  In  cold  weather,  however, 
the  ice  cracks  under  thermal  atnd  wind  stress.  During  storms  ridging 
and  rafting  also- contribute  a  variety  of' noises  like  hammering,  slamming,, 
creaking,  etc.  Blowing  snow  adds  to  the  background  with  a  sound  like 
sandblasting.  Under  these  conditions  the  ambient  level  may  increase  as 
much  as  20  dB. 

Figure  47  summarizes  the  results  of  analysis  of  more  than  100 
ambient  recordings.  A  typical  "quiet"  ambient  (curve  A )  is  .generally 
some  5-10  dB  below  sea  state  zero.  There  is  evidence  of  a  broad  peak 
centered  around  20  Hz.  This  suggests  that  propagation  from-distant 
noisier  areas  is  involved,  since  the  effects  of  scattering  loss  are  quite 
apparent.  Curve  B  is  a  typical  cold  weather  situation  where  the. air 
temperature  has  dropped  suddenly.  With  the  rapid;decline  in  temperature; 
the  ide  comes  under  stress,  and  starts  cracking.  Under  these  conditions, 
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the  noise  background  rises.  Curve  C  is  the  highest  ambient  noise  level 
recorded  during  the  observation  period. 

Ganton  and  Milne16’17  report. ambient  noise  measurements  made 
within  the  Canadian  Archipelago  with  bottomed  hydrophones  in  various 
seasons.  In  the  summer,  the  spectra  fell  within  ~  5  dB  of  the  Knud  sen 
sea  state  zero  curve  above  100  Hz  and  are  "flat"  from  10~100  Hz.  The 
higher  levels  were  associated  with  ice  motion  which  caused  rafting  of  the 
rind  that  sounded  like  "escaping  steam".  In  the  cold  weather  with  de¬ 
clining  air  temperature,  the  cracking  sound  approached  a  steady-state 
hiss  and  increased  the  high  frequency  levels  to  some  20  dB  above  sea 
state  zero.  A  blowing  snow  contribution  to  the  noise  background  during 
windy  periods  was  also  identified. 

Buck?6  also  reported  ambient  measurements  from  20^1000  Hz 
made  in  the  Beaufort  Sea  over  a  two  weeks  period  in  April  1963.  His 
average  value  .far  the  period,  is  about  6  dB  above  sea  .state  zero  with  a 
total  variation  of  roughly  -  10  dB.  The  higher  levels  were  also  associ¬ 
ated,  with  thermal  cra.cking. 

Figure  48  is  a  comparison  of  acoustic  and  seismic  "ambient  noise". 
The  seismic  measurements  were  made  with  a  vertical  component  seis¬ 
mometer  mounted  pn  the  ice  surface^?,  with  the  results  presented  in  terms 
of  a  velocity  frequency  spectrum.  The  acoustic  noise  spectrum  is  derived 
from  the  acoustic  pressure  spectrum  by  u2  =  p2/p0zc02  .  This  implies  that 
the  seismic  disturbances  in  the  ice  and,  probably,  in  the  bottom  also,  are 
coupled  to  the  water  by  direct  acoustic  radiation. 

Occasionally  during; the  summer,  noises  of  obvious  biological  origin 
are  heard.  These  include  long  gliding  whistles,  short  chirps,  braying 
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sounds  and  grunts.  Figure  49  shows  several  Nils silyzer  records  of  these 
noises,  illustrating  the  differing  frequency  characteristics.  An  effort  was 
made  to  identify  these  sounds  with  recordings  made  by  Woods  'Hole  and 
others.  Although  there  was  a  similarity  ,to  the  scmnds  of  certain  whales', 
porpoise,  seals,  etc.  ,  which<are  known  to  frequent  the  Arctic  Ocean,  no 
positive  identification  could;  be  made. 
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FIGURE  7  DISPERSION  TIME  vs  TRAVEL 
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FIGURE  13  MISSILYZER  SPECTRA  (A  SHOWS  BOTTOM  ARRIVALS 
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FIGURE  15  TRANSMISSION  LOSS  vs  RANGE 
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FIGURE  16  TRANSMISSION  LOSS  vs  RANGE 
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FIGURE  18  TRANSMISSON  LOSS  vs  RANGE 
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FIGURE  20  ,  TRANSMISSION  LOSS  vs  RANGE 
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FIGURE  26  TRANSMISSION  LOSS  vs  RANGE 


LEVELS  MEASURED  AT  800  KILOMETERS. 


FREQUENCY  .  HZ 

FIGURE  28  WESTON'S  SOURCE  LEVEL  SPECTRUM 
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FIGURE  32  EIGEN  VALUE  NOMOGRAPH  FOR  BILINEAR,  GRADIENT 
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FIGURE  33  GROUP  VELOCITY  NOMOGRAPH  FOR  BILINEAR  GRADIENT 
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FIGURE  34 
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FIGURE  38  REVERBERATION  GRAZING  ANGLE  vs  TIME 
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FIGURE  44  SCATTERING  STRENGTH  vs  GRAZING  ANGLE 
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FIGURE  46  SYSTEM  NOISE  SPECTRA 
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FIGURE  48  SEISMIC/ACOUSTIC  ICE  SURFACE  VELOCITY  SPECTRUM 
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